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SUMMARY 

The incompressible boundary- layer theory of Truckenbrodt was ex- 
tended using a modified Stewartson-type transformation to include com- 
pressible boundary- layer development over insulated surfaces. The vari- 
ation of turbulent shear stresses with Mach number of the bounding poten- 
tial stream was considered using the reference- temperature method of 
Eckert. The explicit technique of Truckenbrodt for determining the pro- 
file parameter in the incompressible theory was unchanged so that the 
profile parameter as well as the momentum thickness can still be evalu- 
ated by simple quadratures in the present theory. 

The method is applicable to two-dimensional and axi symmetric, lami- 
nar and turbulent boundary layer, and accounts for the pressure gradient 
along the wall for both compressible and incompressible flow. 

Experimental measurements of the boundary layer on small axi sym- 
metric bodies on the wall of a supersonic wind tunnel were compared with 
theoretical predictions. For most cases studied, agreement between ex- 
periment and theory was within 10 percent. 


INTRODUCTION 

It was first shown by von Harman (ref. l) that the differential 
equations describing the motion of the fluid in a boundary layer can be 
simplified by use of integral thickness parameters. The resulting inte- 
gral equations then consider the growth of these parameters (momentum 
and displacement thicknesses) in the external stream direction. Refer- 
ences 2 and 3 further showed that for incompressible flow all velocity 
profiles could be reasonably represented by one parameter family; this 
profile parameter generally being the ratio of displacement to momentum 
thickness. Determination of any two of these three parameters, momentum 
thickness, displacement thickness, and profile parameter, therefore, is 
sufficient to describe the over-all boundary- layer profile characteris- 
tics for most engineering applications. 
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The momentum equation has been solved by the Kamdn-Pohlhausen- 
Holstein method (ref. 2, pp. 93-100) for the laminar incompressible case. 
A method of calculation of the plane turbulent boundary layer was first 
made by Gruschwitz (ref. 4, pp. 93-98). Von Doenhoff and Tetervin 
(ref. 3) also presented a method for calculation of two-dimensional tur- 
bulent boundary layer, the procedure of which was later simplified by 
Garner (ref. 5). Explicit expressions for momentum thickness in terms 
of a simple quadrature were obtained by Buri (ref. 6) and Maskell (ref. 

7) for plane turbulent flow, and by Walz (ref. 8) for plane laminar flow. 

Recent progress has been made by Truckenbrodt (ref. 9) who devised 
a method for computing laminar and turbulent boundary layer for both 
two-dimensional and axisymmetric flow. The main advantage of this method 
over references 2 to 8 is that it enables the computation of the profile 
parameter as well as the momentum thickness in explicit form by use of 
simple quadratures. 

References 2 to 9 all consider incompressible boundary layer. 
Stewartson (ref. 10) presented relations which transformed the two- 
dimensional, compressible laminar boundary- layer equation to the form of 
the incompressible case. Cohen and Reshotko (refs. 11 and 12) applied 
these transformations to the laminar boundary layer including the effects 
of heat transfer. Van Le (ref. 13} suggested that the Stewartson-type 
transformation be used for turbulent as well as laminar flow if time 
average values were taken for the variables. Reshotko and Tucker (ref. 
14) and Mager (ref. 15) used these transformations for analysis of shock- 
induced turbulent boundary- layer separation. 

This paper uses a modified Stewartson transformation to change the 
final equations of Truckenbrodt to compressible coordinates for both 
two-dimensional and axisymmetric flow. The results should then be ap- 
plicable to laminar or turbulent, two-dimensional or axisymmetric, com- 
pressible or incompressible adiabatic flows with or without surface pres- 
sure gradients. The explicit form developed by Truckenbrodt for evalu- 
ating boundary- layer profile parameter and momentum thickness by means 
of simple quadratures is retained. 


SYMBOLS 

A friction term defined by eq. (5) 

a sonic velocity 

b constant in eq. (10) 

constant in eq. (2) 
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local coefficient of friction 

average coefficient of friction 

dissipation (energy converted into heat) 

profile parameter 5*/0 

profile parameter &**/ 9 

parameter used in eq. (9) 

profile parameter (fig. l) 

total length of boundary layer in longitudinal direction 
local Mach, number 

constant a 1 for laminar flow, l/6 for turbulent flow 
pressure 

transformed body radius 
body radius 
temperature 
turbulence energy 

longitudinal velocity component at edge of boundary layer, y s 6 

longitudinal velocity component 

transformed longitudinal coordinate 

longitudinal, coordinate 

transformed normal coordinate 

normal coordinate 

constant in eq. (10) 

ratio of specific heats 

boundary- layer height 
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5* displacement thickness 

5** energy thickness 

e variable in eq. (24) 

£ defined by eq. (23) 

6 momentum thickness 

x defined by eq. (22) 

[i absolute viscosity 

v kinematic viscosity 

£ defined by eq. (l2) 

p density 

t shear stress (includes both laminar and time average turbulent 

stresses if any) 

<p cone half angle 

Subscripts: 

e edge of boundary layer, y « 5 

I incompressible flow in "physical" space 

1 transformed or incompressible flow 

2 laminar 

0 free- stream stagnation conditions 

p flat plate 

r recovery 

s starting, or initial, conditions 

T transition 

t turbulent 

w wall or surface value 
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Superscript: 

based on Eckert's reference temperature used everywhere except 
when modifying 8 


THEORY 

Incompressible Momentum Thickness 

A brief review, to aid in the understanding of the final expressions 
used from reference 9, will be given. 

The Truckenbrodt method of computing momentum thickness is based on 
Wieghardt's kinetic energy equation, 


-JL A. (U^RH. 6. ) m 2 ~ - tJ . 
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instead of the momentum integral equation used by references 2 to 8. Use 
of the energy equation offers the advantage that the term representing 
the sum of the dissipation and turbulence energies (d + t/p Q U^) is al- 
most independent of profile parameter Rj_, whereas the friction coeffi- 
cient, which enters into the momentum equation, varies rapidly with pro- 
file parameter. The dissipation term was studied as a function of Hp 

and U^i/ V q by use of Hartree profiles for laminar flow and the results 
of references 16 and 17 for turbulent flow. 


Equation (l) was therefore integrated, assuming mean values for Hj_, 


resulting in 




A f Vf 
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where 


Aj m 0 \ II 

the integration constant C » R^ +D l — A S J 

The subscript s denotes initial conditions. The constant n 
for laminar flow and l/6 for turbulent flow. 


( 2 ) 


equals 1 


The coefficient A was evaluated by Truckenbrodt by assuming that 
the average shearing stress at the wall was of the same form as that for 
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a flat plate at zero incidence. For a flat plate 
d±; 8 “ 0 if X s * 0, so that equation (2) 1)60011168 



is constant and 


and since for a flat plate from the momentum integral equation 

0 c f,av 

1 ** 2 


-PrlP'H 


(3) 


The average wall friction coefficient was in turn evaluated for laminar 
flow by the Blasius equation 



av 


1.328 



and for turbulent flow by the Falkner expression 



av 




0.0306 



(4) 


(5) 


so that 


A * 0.441 for laminar flow 
A * 0.0076 for turbulent flow 

The two-dimensional case is obtained by omitting R in equations 
(l) and (2). It is interesting to note that the formula for momentum 
thickness (eq. (2)} has the name form as that given by Buri (ref. 6) and 
Maskell (ref. 7) for turbulent flow, and by Walz (ref. 8) for laminar 
flow, in spite of widely different methods of derivation. The exponent 
on Uj0i/vo is given as 0.25, 0.2155, and 0.167 for turbulent flow by 

references 6, 7, and 9, respectively, whereas the exponent on the Uj_ 
term is given as 4.0, 4.2, and 3.33, and the A constants are 0.015, 
0.01173, and 0.0076. For laminar flow, n » 1 in both references 8 and 
9. Reference 8 lists A as 0.470 in place of the 0.441 value in 
reference 9. 
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Wlien considering both laminar and turbulent flow, the final expres- 
sion by Truckenbrodt for momentum thickness may be written as 



where Xp is the longitudinal distance to transition. The expression 

inside the brackets is the contribution of the laminar boundary layer, 
and the remaining part inside the braces is the contribution of the tur- 
bulent boundary layer. 


Incompressible Profile Parameter 

As stated in the INTRODUCTION, the method of reference 9 for deter- 
mining the profile parameter differs from the other methods in that 
Truckenbrodt succeeded in determining the profile parameter in an ex- 
plicit form. The momentum equation 


U£R 


dX^ U i R0 i) + U, dX 


9 1 dU i T w,i 


fo4 


was subtracted from the kinetic energy equation (l) making 


(?) 


dHt 

dT 
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Note that this expression is independent of R. Equation (8) was re- 
arranged to the following form: 
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where the profile parameter 1^ is defined by 
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The quantity Lj_ was arbitrarily chosen as . zero for zero pressure gra- 
dient flows so that the value of Hp^i was set at 2.6 for laminar flow 
and at 1.4 for turbulent flow. 

The data of references 16 to 19 for turbulent flow, and the 
Hartree profiles (numerically evaluated in ref. 20) for laminar flow were 
used to arrive at approximate expressions for K as a function of L-^, 
and Lj^ as a function of Ej_. 

An adequate approximation for the quantity K(L^) is 

K(L ± ) » a(Li - b) (10) 


where 


a 


6 

6 . 


87 for laminar flow with pressure drop 
53 for laminar flow with pressure rise 

0304 for turbulent flow 

for laminar flow 

0305 In U^i/vQ - 0.23 for turbulent flow 


The resulting relations between Lj_ and are shown in figure 1. 

Substitution of equation (10) in equation (9) results in a first 
degree linear differential equation, then solved for L^. The result is 
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for laminar flow 
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for turbulent flow. The parameter Is defined by the two equations 

(l2a) and (12b) since the change of profile parameter through the tran- 
sition region will be handled as a discontinuity, as will be discussed 
later in the Initial Conditions section. 

The calculation procedure may be briefly summarized as follows: for 

a given velocity distribution over a known body, the momentum thickness 
distribution is first calculated using equation (6); then the profile 
parameter Lj_ is calculated as a function of the longitudinal distance 
from equations (ll) and (12), and finally, the conventional form factor 
Hi is obtained from the values of using figure 1. 

Transformation to Compressible Case 

A modified Stewartson transformation for a Prandtl number of one 
and an adiabatic flow (see Appendix and symbol list) is 



Where x, y, and r are the compressible "physical" coordinates and 
X, Y, end R are the transformed coordinates. Use of these transforma- 
tions in either the Karman momentum equation or the Wieghardt kinetic 
energy equation yields the following relation for the friction coeffi- 
cient (see Appendix) . 
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Therefore the transformations satisfy the known relations for the lami- 
nar friction coefficient. 


For turbulent flow* the total shear term includes turbulent as well 
as laminar shear stresses , and no simple expression for the total stresses 
is available. Eckert (ref. 21), however, found that the form of the ex- 
pression for a flat plate friction coefficient can be made invariant with 
the Mach number of the external stream when the properties of the fluid 
are based on a proper reference temperature. This reference temperature 
is described as 

T*« T e + 0 . 5 (T w - T e ) + 0.22(T r - T e ) 
or 

m* T w 

i— => 0.72 — + 0.28 

T e T e 



Expressions for friction coefficient obtained for incompressible flows 
can then be extended to apply to compressible flows. From reference 21 


c f 




( 15 ) 


Applying this method to equation (4) of this report and assuming that 
the viscosity varies linearly with temperature, 
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The effect of compressibility on the turbulent friction coefficient could 
be shown through the ratio 



where c.p 
r I 


is defined as the constant fluid- property (incompressible) 


friction coefficient for the same Reynolds number as the compressible 
stream and should be kept distinct from the incompressible friction co- 
efficient in the transformed plane. The agreement between the re- 

sults of equation (16) (y = 1.4) and those of references 21 to 26 is shown 
in figure 2. The curve of equation (16) is somewhat lower than that of 
reference 21, based upon the Schultz -Grunow equation (ref. 27). Closer 
agreement of equation (16) with reference 21 and the experimental data 
could have been obtained by adjustment of the constants in the expression 
for reference temperature. This adjustment was not believed necessary 
for purposes of the present study. 


A transformation of the compressible friction coefficient by use of 
equations (Al) , (5), and (15) results in 


C ft 
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c f 



0.0306 p* 




( 17 ) 
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In figure 3 it is shown that 


a * \®0 / 


6/7 


so that equation (17) becomes 


c f + 


(•o) 


88/49 


0.0306 


m 


(18) 


Therefore, the transformed friction coefficient may be expressed as 



Recall that friction coefficient enters into equation (6) through the A 
terms, where for turbulent flow 



A combination of this equation with equation (19) yields 

/vf / 21 

At = 0.00761 

When Mq -*■ 0, — -*• 1 so that equation (6) remains unaltered for incom- 


( 20 ) 


pressible flow. 


Therefore, under the circumstances considered (insulated surface 
and Prandtl number equals one) the laminar boundary layer is exactly 
transformed to the incompressible form, whereas the turbulent boundary 
layer is transformed to a form which reduces to the incompressible form 
for a Mach number of zero . 
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A transformation of equation (6) by means of equations (Al) and (20) 
to get the compressible flow momentum thickness gives 
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( 21 ) 


( 22 ) 


(23) 


For the completely laminar boundary- layer case, x equals £ and 

for the completely turbulent case, x equals only the second term on the 
right side of equation (22) . 


When determin i ng the profile parameter, it is convenient first to 
determine Iq and then find Hj_ and H. By substituting equations (Al) 

and (20) in expressions (ll) and (12) , and by using the definitions of 
a and b, the following equation results: 
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where 

6 5 

6a ^ for laminar flow with pressure drop 

g 

e as £ for laminar flow with pressure rise 


6 



for turbulent flow 


b a 0 for laminar flow 


b = 



0 . 23 for turbulent flow 


Figure 1 can now be used, to determine A combination of equations 

(A3) and (A5) results in the following expression for H: 


H = 



(25) 


Boundary- layer profile and momentum thicioiess for compressible flow 
can thus be determined in explicit form by use of equations (2l) to (25). 
Recall that these equations are written in the more general or axisym- 
metric case. For two-dimensional flow merely omit r whenever it 
appears . 


Initial Conditions 


Initial values for laminar flow are derived in reference 9 as 
follows: For two-dimensional stagnation- point flow. 


^i,s “ 0 


*i,s 


.0260 


0.271 




or 


0.271 ^ a e~y 

it 1 duW)/ 


9r-i3 

2r-2 
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-i 


For axisymmetric stagnation-point fl ow, 

Lj^ s = 0.0195 
0.235 


9 ±,b ~ 


’ a[TeL 

Vo ax 


or 

9r-i5 
2 


For both, two-dimensional and axisymmetric flow over a sharp-edged or 
pointed body., 


0.235 



a - 0 


and 


(26) 


e i,s “ 0 s “ 0 

There will be a relatively small change in 0 through the transi- 
tion region if the length of the boundary- layer travel is large, there- 
fore, it is reasonable to assume that 9 ^ at the end of the transition 

is equal to 0, at the start of transition. Appreciable change in the 
profile parameter, however, may take place through the transition region. 
For incompressible flow with no pressure gradient, reference 9 lists the 
measured results of references 27, 16, and 17 as shown in figure 4. These 
results show only a small variation of change of the incompressible pro- 
file parameter through transition as a function of the Reynolds number. 
Some data for transition showing a large effect of the pressure gradient 
may be found in reference 7 for incompressible flow. Inaccurate assump- 
tions of the initial value of H are quite permissible, however, be- 
cause the influence of the values selected for an initial condition van- 
ishes rapidly with the distance downstream of the transition. This will 
be demonstrated later in the EXPERIMENT section. 


EXPERIMENT 

The method advanced in this report was compared with a limited 
amo unt of experimental measurements of the boundary layer on several 
small axisymmetric bodies and on the walls of a large supersonic tunnel. 
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The measurements were made with small pitot tubes. Momentum thickness 
and form factors were computed with the assumption that static pressure 
and total temperature were constant in the y-direction through the 
boundary layer. These assumptions are the same as those of the appendix. 
The Reynolds number (UZ/v) was approximately 5x10® for the axisymmetric 
bodies and approximately 2x10® for the supersonic tunnel. The transition 
Reynolds number was approximately 1.5x10® on the axisymmetric bodies. 

The boundary layer on the tunnel was essentially completely turbulent. 


Comparison of theory with measurement made on cones having half 
angles of 20° and 25° is shown in figure 5. For the case of the constant 
pressure gradient (dp/dx w 0), equations (21) can be integrated to give 



Thus, no integrals need be evaluated to determine the momentum thickness. 
The largest difference between experimental and theoretical momentum 
thickness is 7 percent, which is probably close to experimental accuracy, 


The change of the profile parameter (H or H-j_) with the turbulent 
boundary- layer growth is negligibly small on these cones since there is 
no pressure gradient in the x- direction and the length of the boundary- 
layer development is approximately 1 foot. These data therefore were used 
as a check of the profile-parameter transformation (eq. (25)). The com- 
pressible profile parameter was computed from the incompressible profile- 
parameter values of 1.3 and 1.4. This is the maximum range of % ob- 
tained using equation (26) and figure 4. Good agreement with the ex- 
perimental trend was obtained. 


Figures 6 and 7 compare theoretical and experimental results for 
the cases of moderate and high, adverse pressure gradient (dp/dx > 0) . 
The abrupt pressure gradient on figure 7 is due to an externally gener- 
ated shock striking the body. These bodies are all conical for 
0< 0.233, in which range it was assumed that dp/dx = 0. The ex- 

perimental Mach number variation with distance in the x- direction is 
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shown in figures 6(a) and 7(a) , and the calculated momentum and dis- 
placement variation is shown in figures 6(b) and 7(b). The agreement is 
f generally within 10 percent, which is again close to the accuracy of the 

data. 

The effect of making an erroneous assumption of the initial value 
of the profile parameter is also shown in figure 6(b) by the line made 
up of short dashes. An initial value of 2.2 (near separation value) was 
assigned to even though the pressure gradient in this region was 

small. In spite of this poor assumption, the solution converged to the 
original curve (in which was assumed to be 1.4) in a very short 

t distance. 

s 

The case of a favorable pressure gradient (dp/dx < 0) on the wall 
of a supersonic tunnel is shown in figure 8. Since the flow in this 
case was neither completely axisymmetrie or two-dimensional, an effec- 
tive radius equal to the tunnel wetted perimeter at each axial station 
and divided by 2it was used for r in equation (2l) . This computed 
result was compared with the experimental 0 and H measured on two ad- 
jacent tunnel walls. Experimental agreement with theory was generally 
within 10 percent. 

Only small changes in the incompressible profile parameter (%) were 
. obtained in the cases of the favorable and moderate adverse pressure 

gradient. Large changes, however, were obtained for the compressible 
profile parameter (H) . Therefore, for rough approximations in these 
cases and when initial conditions are known quite well, it may be per- 
missible to hold Hi constant and determine H from equation (25) only, 
thus eliminating the use of equation (24) . Obviously, the elimination 
of equation (24) would not be permissible for boundary layers in the 
region of large adverse pressure gradients. 


CONCLUDING REMARKS 

A method has been presented herein which enables a determination of 
the compressible boundary- layer growth and profile over both two- 
dimensional and axisymmetrie bodies. Both momentum thickness and profile 
parameter are obtained by simple quadratures. The method includes both 
laminar and turbulent flow; however, the experimental check involved 
mainly turbulent flow. The momentum thickness at the transition point 
was calculated to be approximately 10 percent of the total momentum 
thickness for the cones and axisymmetrie bodies. Comparison was made 
between theory and experiment for flows with zero, adverse, and favorable 
pressure gradients in the direction of boundary- layer travel. For most 
cases studied, agreement between experiment and theory was within 10 
percent . 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, April 16, 1957 



18 


MCA TN 4022 


APPENDIX - TRANSFORMATION OF MOMENTUM AND ENERGY EQUATIONS 


The modified Stewartson transformation used in this report is: 

„ 3]tl ^ 

r-i 


x = 




dx 


Y = — J5_ dy 

a °Jo p ° 


u i - a 0^e 


ao 


U 


> (Al) 


Uj = U 

1 a e 

and for axisymmetric flow it will he shown that the radius is 
invariant with these transformations so that 

R *» r 


J 


These relations will he used to transform the compressible momentum 
and energy equations to corresponding equations of the incompressihle 
form. In like manner , it could he shown that these relations can he used 
for the Inverse transformations, that is, from the Incompressible to the 
compressible form. 


It will he assumed that there is no heat transfer, that the Prandtl 
number equals one, that static pressure is constant in the y-direction 
for 0 4 y^: 6, and that the flow external to the boundary layer Is 
Isentropic. 


Any of the equations in this appendix can he applied to two- 
dimensional flow by omitting the r and R terms. 




Momentum Equation 

The Karman momentum integral equation for axisymmetric compressible 
flow can he written as follows: 


1 _d_ 

rp e U 2 **■ 


<***> + S f 


6* 



(A2) 


4326 



CH-3 back 4326 


NACA TN 4022 


19 


By 


To transform this equation* some useful relations will "be derived, 
definition 



and by using equation (Al) 


then also 


By definition 



(A3) 


with equation (Al) 



r+i 

(AS) 
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From equation (Al) 


and 


dU. 

dX 


4r-2 

r-i 



Hi) s*' 

_ T - l faAMl dU 
dx 2 \a e j U dx 


5y-l 
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since 


r-i 


^o 

8,0 



(A6) 


and 


dp e - -p e U dU 


(A7) 


Therefore , 


4r-2 

faA r " 1 dUi 
dU V^i dX 

i + T ~ - m| 


(AS) 


and 


<a 


dX 


y - 1 M e a e ^i 

2 1^ diT 

i -r g 


(A9) 


By substituting equations (Al) , (A3)> (A5), (A6), and (AS) in equation 
(A2) 
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80 d ^ ae> )Ui R0 ii + 
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After using equation (A9) and collecting terms, there results 

2 nr- 


4 + £ ®r - ($)' 


w 


Pe u " 


(A10) 


and thus the momentum equation takes the incompressible form if 
(ao/a e ) 2 T w /p e U 2 = ^± }V /pqP±’ This transformation of the friction term 

is discussed in the THEORY section of this report. 


Energy Equation 


The kinetic energy integral equation of Wieghardt for axisymmetric 
compressible flow may be written in the following integral form: 


,8 


-8 




2p_U°r 


dy 


(All) 


Equations (Al), (A6) , and (A8) are substituted in equation (All) to give: 



p O a O 

o,- 

p 1 

^e e 
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By using equations (A4) and (A9) in equation (A12), 


2U^R dX 




1 M e 


2 U A dX 

nxzr^i 






and by collecting terms , 


1 _d_ 

2ufl? ax 


(U?R5f) 



(A13) 


which is the incompressible form of the Wieghardt kinetic energy equation 

if (ag/a e ) 2 Tr w /p e U 2 = ^^ w /pqU? . Thus, if the friction transformation is 

satisfied in the momentum equation, it is satisfied in the Wieghardt 
kinetic energy equation. 
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Mach number, M e 

Figure 2. - Effect of compressibility on skin friction of a 
flat plate. 


4326 



Mach number. Mg 


Figure 3. - Comparison of speed of sound based on reference 
temperature of Eckert with speed of sound based on stagna 
tion temperature. Ratio of specific heat, 1.4. 
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Figure 4. - Change of incompressible form factor in the range of transition from laminar to 
turbulent flow (ref. 9). So pressure gradient. 
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(a) Mach number gradient. 


Figure 6. - Boundary layer on an axi symmetric body with mod- 
erate adverse pressure gradient. 
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Figure 8. - Boundary layer on the vails of a wind tunnel vith favorable pres- 
sure gradient. 
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Figure 8. - Concluded. Boundary layer on the walls of a wind tunnel with favorable urea 
sure gradient. * 
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